One-qubit reduced states of a pure many-qubit state: polygon inequalities 
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We show that a necessary and sufficient condition for a set of n one-qubit mixed states to be 
the reduced states of a pure n-qubit state is that their smaller eigenvalues should satisfy polygon 
inequalities: each of them must be no greater than the sum of the others. 
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In this paper we study the quantum analogue of the 
marginal distributions of a joint probability distribution. 
The results reveal some surprising aspects of pure states 
of multipartite quantum systems. 

In a pure state of many particles, each subset P of the 
particles is in a mixed state /Op, the reduced state of the 
subset P. These reduced states are subject to conditions 
such as the consistency conditions 



trg Ppuq — tr^ ppuR 



(1) 



for any subsets P, Q, R such that P and Q are disjoint 
and P and R are disjoint. However, equations of this 
type do not exhaust the conditions: Linden, Popescu and 
Wootters have given an example of a set of three two- 
qubit density matrices which satisfy all such conditions 
but cannot be the reduced two-qubit states of a pure 
state of three qubits. Another example is given by the 
system of four qubits, in which it is impossible for all 
three reduced two-qubit states to be totally mixed 
though all the consistency conditions are satisfied. 
The full set of conditions is not known. 

In this paper we completely determine the possible 
one-qubit reduced states of a pure state of n qubits. 
(There is no problem for mixed states, since any n one- 
qubit mixed states pi, . . . pn are the reduced states of 
the n-qubit mixed state pi O • ■ • ® Pn-) This is equiv- 
alent to determining the set of possible eigenvalues for 
each qubit, i.e. the marginal probability distributions. 
Brun and Cohen Q have pointed out that in the case 
n — i the determinants of the reduced states must sat- 
isfy triangle inequalities, but not every set of three 2x2 
density matrices satisfying these inequalities is a possible 
set of reduced states of a pure three-qubit state. We will 
show that such density matrices must satisfy the stronger 
condition that their smaller eigenvalues satisfy triangle 
inequalities. For n qubits, a necessary condition is the 
obvious generalisation of the triangle inequalities: 



Afe < Ai + ■ • • + Afe_i + Afe+i + • • • + A„ 



(2) 



where the eigenvalues of the density matrix pk of the fcth 
qubit are Afc and 1 — Afc with A/t < 1/2. These inequalities 
completely characterise the possible sets of eigenvalues, 
and therefore the possible reduced one-qubit states, of a 
pure n-qubit state. 



These results reveal some surprising aspects of pure 
quantum states. One concerns the comparison with the 
classical notion of a pure state in the sense of a prob- 
ability distribution. The notion of a pure state applies 
to both quantum and classical systems, if we consider a 
state of a classical system to be a probability distribution 
for the values of variables which actually have precise val- 
ues. Then a pure state of a classical system (defined to 
be an extreme point of the convex set of probability dis- 
tributions) corresponds to perfect knowledge of the vari- 
ables, so that all probabilities are or 1. A pure state of 
a quantum system corresponds to maximum knowledge, 
but the characteristic feature of quantum mechanics is 
that even in a pure state there are variables for which 
perfect knowledge is not available; indeed, any possi- 
ble probability distribution will occur for some variable. 
However, given any pure state there is always some phys- 
ical variable for which the probabilities are all or 1 (for 
example, the projection operator onto the state). This 
is the quantum version of the above characterisation of 
classical pure states. 

This analogy between quantum and classical pure 
states disappears when we consider marginal states. 
Classically, a joint distribution of many variables is pure 
if and only if each one- variable marginal distribution is 
pure. We know that this is not true in quantum me- 
chanics; the fact that a pure multipartite state may give 
mixed one-party states is the phenomenon of entangle- 
ment. Nevertheless, our results place limits on the one- 
party states; not every probability distribution is possi- 
ble. 

The most surprising aspect of this is that it is a state- 
ment about an overall property of the whole system, 
namely the purity of its state, which depends only on 
local measurements. We emphasize that this is only a 
negative statement: from the knowledge that the one- 
party reduced states satisfy the polygon inequalities one 
cannot normally deduce that the multipartite state must 
be pure (this is only possible if the one-party states are 
themselves pure); but the local information in a viola- 
tion of the polynomial inequalities does give the overall 
information that the multipartite state cannot be pure. 

To prove the assertion (O , we expand the pure n-qubit 
state ll-) in the Schmidt basis Let l^*"^) and 
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b^i^) be the eigenstates of the reduced state pk with 



eigenvalues Xk and 1 — Xk respectively; we will show that |$]^ ) in Q), 



Ai < A2 + --- + A„. 



(3) 



We assume Ai > (if not, the inequality ||2Jl is trivially 

satisfied). We note that I^q"'^'') and \(t>'"i^) are the states of 
the first qubit occurring in the Schmidt decomposition 
of l^*) when the system is divided into two parts, one 
containing the first qubit and the other containing the 
rest, so 

|*)=^|0«)|$W)+i?|0i^))|$W> (4) 

where I'&g^'') and l^^^"*) are orthogonal normalised (n— 1)- 
qubit states, A and B are real and — Xi = 1 — . 
Since Ai is the smaller eigenvalue, < A < B. We write 
\ii . . . in) = ■ • • and expand l^*) in this basis: 



m^A 



E 



B 



with 



I2,l3,...,ln I2,l3,---,l 



Jl«2«3 ■■■'in) 



E \b^.^3■■■^J' = ^ ■ (5) 



Then 



Xk = A^ E 



i2,i3,---,in;ik=0 

12,13, ■■■,in;ik=0 

SO that 

X2 + ■ ■ ■ + Xn = A^ -^i2i3---in ki2i3---in 1^ 



B' J2 N,,,,...,Jh, 



l2«3---«n I J 



where 



Ni2i3-t„ = n - 1 - (Z2 + 13 H h in) (6) 



is the number of times "0" appears in (^2, ^3, ■ • • , in)- 
Hence Ni^i^...,^^ > 1 if («2, «3, ■ • • , ««) 7^ (1, 1, • ■ ■ , !)■ This 
implies that 



Xn > A'^ \ ^ |a»2»3- - 




= ^^(l-|aii...ir) + B"(l-|6ii...i 



> A^{2-\an-i\'-\bii...i\^), (7) 



since B^ > A^. Now by the orthogonality of l^'^^) and 



^ ^ ^i2i3---in bi2i3---in — 0, (8) 
i2,i3,---,in 

so, using the Cauchy-Schwarz inequality, 



|an...in6ii...iP < E 



1 l2,«3,---,Jn 



E l^«2»3-i„P - l^ll-l| 

= (l-|an...i|')(l-|6ii...in. (9) 

Therefore 

|aii...iP + |6ii...iP < 1. (10) 

Hence from Eq. Q we have 

A2 + A3 + ■ • • + A„ > ^2 ^ Ai . (11) 

Clearly there are similar inequalities with each of the 
Xk on the right-hand side. We call these the polygon 
inequalities. 

To show that these inequalities define exactly the set 
of possible one-qubit reduced states, we prove that given 
any real numbers {Ai, . . . , A„} lying between and 1/2 
and satisfying jSJ we can find an n-qubit state for which 
Xk is the smaller eigenvalue of the reduced density matrix 
Pk- We will suppose that Ai is the largest of the numbers, 
so that 

0<Ai<Ai<l/2, i = 2,3,...n. (12) 

For n = 3 let 

1^-3) = a|100) + 6|010) + c|001) + d|lll) , (13) 

where a, b, c and d are real and satisfy a'^+b'^+c^+cP = 1. 
Then for each of the three qubits the eigenvectors of the 
one-qubit reduced state are |0) and |1). Let A^ be the 
eigenvalue corresponding to |0) for the zth qubit. Then 



Ai = b^ + c^, 
A2 = c^ + a^, 
A3 = a2 + 62. 



These can be solved as 



= ^(Aa + Ag-Ai) , 
6^ = 1 (A3 + Ai - A2) , 
= i (Ai + A2 - A3) . 



(14) 
(15) 
(16) 



(17) 
(18) 
(19) 
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Thus if the Xi satisfy the triangle inequahties there is a 
real solution {a,b,c) of the equations, and if Xi < 1/2 
this solution satisfies 

a2 + 52 + = ^^ + ^^ + ^^ < 3/4 , (20) 

so that + + < 1 and it is possible to find a nor- 
malised state 1^1'). Hence, there is a state with arbitrary 
eigenvalues satisfying Eqs. ((2Jl and (|12|) for n — 3. 

Now assume that this is true for any (n — l)-qubit 
system. Suppose that A„ is the smallest of the A^, and 
define Ai = Ai — A„. There are two cases to consider. 
If Ai > A,; for 2 < i < n - 1, then (Ai, A2, . . . , A„_i) 
satisfies Eqs. Q and H12|l with n replaced by n — 1 and 
Ai by Ai. In the second case, Ai < Xm for some m, say 
TO = 2; we can suppose that A2 > A; for i — 3, . . . ,n — I. 
Then 

Ai + A3 + • • • + A„_i = Ai + (A3 - A„) + A4 + • • • + A„_i 
> Ai > A2 . 

Hence the set (A2, Ai, A3, ... , A„_i) satisfies Eqs. (O and 
(|12|l with n replaced by n — 1 and Ai and A2 replaced by 
A2 and Ai, respectively. In either case there is a state for 
which Ai, A2, . . ., A„_i are the smaller eigenvalues of the 
reduced density matrices. Let this state be 

l*n-l> = |0)|(/') + |l)|^), (21) 

where \(f>) and are (n — 2)-qubit states satisfying 
(0 I 0) = Ai, (1/' I V') = 1 - Ai and {(|)\^p) = 0, and A2, 
A3, . . ., A„_i are the smaller eigenvalues of the one-qubit 
reduced density matrices of \(p) {4>\ + lip) ('01 • Now consider 
the following n-qubit state: 

= |O)|0)|l)+sinx|O)|7/;)|O)+cosx|l)|0)|l). (22) 

The smaller eigenvalue A^ of this state for the ith qubit 
with 2<z<n— lis again the smaller eigenvalue of the 
ith one-qubit reduced density matrix of |(/')(0| -I- I?/') ('01. 
Hence A; = Ai for 2 < i < n — 1. The nth and 1st 
eigenvalues corresponding to the eigenvectors |0) are 

A„ = sin2x(0|0), (23) 
Ai = (0|0)+sin2x(0|0) = Ai + A„. (24) 

Since > 1/2, one can choose x so that A„ ~ A„ . 

Then gives Ai = Ai. 



We conclude by induction that a state with the re- 
quired one-qubit eigenvalues exists for any n. 

Of course, it does not follow that if an n-qubit state 
has one-qubit reduced states whose eigenvalues satisfy 
the polygon inequalities, then the n-qubit state must be 
pure. In general, a given set of one-qubit reduced states 
can be obtained from many different n-qubit states, most 
of which will be mixed: for example, the tensor product 
of the one-qubit states, which will be mixed if any of the 
one-qubit states are mixed. It is only when all the one- 
qubit states are pure that the n-qubit state from which 
they arise must be pure. 

Purity of the n-qubit state places no restriction on 
the eigenstates of the one-qubit reduced states, since any 
pair of orthogonal eigenstates can be transformed to any 
other by local unitary operations. The complete set of 
one-qubit reduced states which can arise from a pure n- 
qubit state is therefore determined by the set of possible 
eigenvalues that we have described. The situation for the 
reduced states of larger subsets is likely to be more com- 
plicated. The example of two-qubit subsets of a system 
of four qubits shows that polygon inequalities ^ for any 
combination of eigenvalues are not sufficient, for these 
would allow all reduced states to be totally mixed (all 
eigenvalues equal to 1/4), which is not possible 2]. It 
would be interesting to know the exact set of marginal 
two-qubit probabilities for this system. 

The situation for qudits (particles whose state spaces 
have dimension d) also appears to be more complicated. 
An argument similar to the above shows that the eigen- 
values of the one-qudit reduced states must satisfy Q 
with each Ai replaced by the sum of all but the largest 
eigenvalue of qudit i (so the largest eigenvalues satisfy 
the same inequality as for qubits), but there appear to 
be further inequalities that must be satisfied. This is 
under investigation. 
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